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From a kinematical point of view, the geometrical informa- 
tion of hamiltonian chaos is given by the (un)stable directions, 
while the dynamical information is given by the Lyapunov 
exponents. The finite time Lyapunov exponents are of par- 
ticular importance in physics. The spatial variations of the 
finite time Lyapunov exponent and its associated (un)stable 
direction are related. Both of them are found to be deter- 
mined by a new hamiltonian of same number of degrees of 
freedom as the original one. This new hamiltonian defines 
a flow field with characteristically chaotic trajectories. The 
direction and the magnitude of the phase flow field give the 
(un)stable direction and the finite time Lyapunov exponent 
of the original hamiltonian. Our analysis was based on a 1^ 
degree of freedom hamiltonian system. 

PACS numbers: 05.45. -fb 



I. INTRODUCTION 

Many systems that are chaotic such as the stochastic 
magnetic field line in toroidal plasma confinement de- 
vices have a hamiltonian representation . The KAM 
theorem [|| deals with the onset of chaos in hamiltonian 
systems, and the description of the kinematics of devel- 
oped chaos involves the ergodic theorem Q . Hamiltonian 
dynamics naturally give rise to a differentiable dynami- 
cal system where the multiplicative ergodic theorem of 
Oseledec is apphcable Oseledec's theorem gives two 
aspects of chaos. First the sensitivity of the dependence 
on initial conditions is measured by various infinite time 
Lyapunov exponents, i.e., the dynamical information of 
chaos. Secondly the characteristic directions associated 
with these Lyapunov exponents give the geometrical as- 
pect of chaos. That is: if two points are separated along 
different directions at the initial time, they can diverge 
or converge exponentially at different characteristic Lya- 
punov exponents. 

Ergodic theory treats the time asymptotic limit, in 
which the infinite time Lyapunov exponents are constants 
and the characteristic directions are functions of position 
only. For finite time, there is a convergence issue for 
both the Lyapunov exponents and their associated char- 
acteristic directions. The convergence of the characteris- 
tic directions is exponential, so the geometrical aspect of 



chaos at finite time is well described by its time asymp- 
totic limit. This is not the case for Lyapunov exponents. 
Finite time Lyapunov exponents suffer from a notoriously 
slow convergence problem Its spatial and temporal 
dependence was discussed by Tang and Boozer who 
gave a direct link between the convergence function in 
finite time Lyapunov exponent and the geometry of the 
vector field defined by the corresponding characteristic 
direction. It should be emphasized that most applica- 
tions of practical interest involve a finite duration of time. 
Hence the finite time properties of chaos, rather than the 
asymptotic properties of chaos, are of real concern. 

In |7| we found that the spatial variation of the finite 
time Lyapunov exponent and the corresponding charac- 
teristic direction are not independent of each other. The 
exact relation was shown in and is restated later in 
this paper, equations (0,^,^). Besides its practical im- 
portance such as those in transport studies 0, it also 
leads to a conceptual advance in the understanding of 
chaos. The finite time Lyapunov exponent (e.g. A) and 
its characteristic direction (e.g. Sao) can be described by 
another hamiltonian in the same phase space with the 
same number of degrees of freedom as the original hamil- 
tonian. If one constructs the corresponding vector field 
(necessarily divergence-free) associated with the phase 
space trajectory of this new hamiltonian, the magnitude 
of this vector field gives the local finite time Lyapunov 
exponent while the direction of the vector field gives the 
corresponding characteristic direction. The new hamil- 
tonian is also chaotic in the same region as the parent 
hamiltonian. One could characterize the chaos in this 
new hamiltonian by invoking a third generation hamil- 
tonian. Consequently, a hierarchy of hamiltonian can be 
bootstrapped from the original hamiltonian of a conser- 
vative system that is chaotic. 

II. GLOBALLY DIVERGENCE-FREE FIELDS 
AND HAMILTONIAN MECHANICS 

A large class of conservative systems which exhibit 
chaotic behavior has a hamiltonian representation. Two 
of the well known examples are the magnetic field B and 
the velocity field v of a divergence-free fluid |^ . An arbi- 
trary divergence-free vector G(a;, y, z) can be written in 
the so-called canonical form. 
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G = V-0 X V6' + V0 X Vx- 



(1) 



The function x('0, 6*, </>) is the hamiltonian of the G field 
fines. Tfie G field line is the trajectory x(t) given by 
equation 



(2) 



To find the G field line in i/', 6, (p coordinates, one has to 
invert the transformation equation x{ip,9,(l)). Combin- 
ing the transformation equation x(?/', 0, 4>) and the field 
line trajectory, equation (H), one arrives at the familiar 
Hamilton's equations. 
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where x is the hamiltonian, 9 the canonical position, 
ip the canonical momentum, and (p the canonical time 
1^. Since continuous transformation, such as x{ip,9,(p), 
preserve topological properties, questions regarding the 
integrability of the field line are then answered by the 
hamiltonian x(V'j 4') alone. 

It should be noted that only a globely divergence-free 
field can be represented by the form of equation (|l|), a 
requirement set by Poincarc's lemma Furthermore, 
the divergence- free field G should not vanish in the region 
of interest. 

A hamiltonian representation of such a vector field is 
desirable since an array of well developed techniques from 
hamiltonian mechanics are available. For example, the 
hamiltonian description of a magnetic field has played a 
major role in the theory of toroidal plasmas 

Before we proceed to the next section which is on the 
ergodic theorem, a digression on the applicability of the 
multiplicative ergodic theorem to hamiltonian chaos is 
useful. A hamiltonian flow, even when it is chaotic, can 
not support a nonuniform invariant density on its ergodic 
subcomponent in phase space. However, numerical sim- 
ulations appear to offer a contradictive picture: there are 
concentrations, or spikes, of the phase space density near 
the boundary of chaotic zones One might question 
the ergodicity of the chaotic sea (the irregular component 
[^) and hence the applicability of Oseledec's theorem. A 
recent study by Meiss in the case of the standard map 
appears to reconcile this conflict. Meiss found that 
these density spikes are transient and they can be ex- 
plained by a Markov tree model taking into account the 
effects of islands-around-islands structure. Hence earlier 



numerical findings |10| are not in contradiction to an in- 
variant measure, and the phase space density can be uni- 
form on an irregular component of hamiltonian systems. 



^In this extended 'phase' space, there are three independent 
coordinates and hence three Lyapunov exponents. 



Despite the lack of a rigorous theory g], these new nu- 
merical evidence and the simple Markov transport model 
put forwarded by Meiss should boost our confidence in 
the ergodicity of hamiltonian chaos. 



III. MULTIPLICATIVE ERGODIC THEOREM IN 
LAGRANGIAN COORDINATES 

The multiplicative ergodic theorem of Oseledec ||] 
complements KAM theorem Q in understanding hamil- 
tonian chaos, especially the limiting case in which the 
regular components occupy a small portion of the phase 
space. This theorem can be understood in terms of a 
general coordinate transformation pl| between ordinary 
space x and the Lagrangian coordinates ^, a widely used 
coordinate system in fluid mechanics. A point with La- 
grangian coordinates ^ is related to x in ordinary space 
by the integral curve of dx./dt = G(x) from initial time 
to to later time t with the initial condition x(^, to) = ^. 
For a pair of flxed (to, t), there is a one-to-one mapping 
between ^ and x(^,t). There is freedom in the choice of 
Lagrangian coordinates ^ due to the arbitrariness of to, 
a point that we will revisit in section [V. Of course, 
to should be fixed for any chosen set of Lagrangian co- 
ordinates, to is usually set to zero for convenience of 
bookkeeping. 

In addition to a set of coordinates, one needs a met- 



ric to specify the physical distance between two neigh- 
boring points. Let x be ordinary Cartesian coordinates, 
then its metric is a unit matrix, dP — dx ■ dx ^ dx^dx^ . 
The same differential distance can also be specified in 
Lagrangian coordinates, dP — gijdS^'^d^^ , where gij = 
(9x/9^* • dx/d^^ is the metric tensor of the Lagrangian 
coordinates ^. Since x(^, 0) = ^, one could also under- 
stand {dxY — gijdCd^^ as an equation which relates the 
initial separation d^^d^^ between two neighboring points 
to their later separation (dx)^. Once this is understood, 
the interpretation of the metric tensor gij in Oseledec's 
multiplicative ergodic theorem becomes transparent, i.e., 
gij is the Oseledec matrix Ay (for its definition, see 
ori). 

Instead of taking the limit of lim(^oo(Aij)^/^* = A and 
diagonalizing the matrix A, we diagonalize the matrix gij 
first and then take the limit of t — > oo. Since gij is a 
positive definite and symmetric matrix, it can be diago- 
nalized with positive eigenvalues and real eigenvectors, 

gij = Aiee + A,„mm -I- A^ss, 

with A/ > A„i > As > 0. There are three Lyapunov 
characteristic exponents associated with vector field G, 
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For a divergence-free field with a hamiltonian represen- 



tation, A? 



-Af > and A° 



0. These are usu- 



ally called infinite time Lyapunov exponents and they are 
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constants in one ergodic region. The eigenvectors e, ih, s 
have weh-defined time-asymptotic hmits as well, but they 
are position dependent, eoo(^), moo(Oj §oo(C)- There is a 
stable manifold theorem in the vicinity of an arbitrary 
point ^ 1^. The §00 (0 is the tangent vector of the local 
stable manifold associated with the Lyapunov character- 
istic exponent A^, the so-called stable direction. There 
is also an unstable direction which is tangent to the lo- 
cal unstable manifold. The unstable direction is defined 
using the time-reversed dynamics. 



in which = (dx)2 = glj^'d^'d^i and glj' = 

A|"^e(-)e(-) + AL"'m(-)m(-) + A^^s(-)s(-). The lo- 
cal stable manifold for the backward time dynamics de- 
fines the local unstable manifold for the forward time dy- 
namics. Two linearly independent vectors define a plane. 
The intersection of two planes gives rise to another vec- 
tor. One can use the intersection of plane (§00, liioo) and 
plane (sio , liiL ' ) to construct another vector m. For the 
case we are interested in this paper, = and m is 
tangent to the center manifold. In the general case 
could be nonzero, m is then the tangent vector of another 
local (un)stable manifold, corresponding to a (positive) 

negative A^. In either case, these three vectors Soo,s^\ 
and m span the three dimensional space. 

The negative Lyapunov exponent means neighboring 
points converge exponentially in time, while a positive 
Lyapunov exponent means neighboring points diverge ex- 
ponentially in time. The ones of most importance are 
the most negative and the most positive Lyapunov expo- 
nents. In applications, one is concerned with the finite 
time Lyapunov exponents. In the case studied in this 
paper, there are two of them and they are defined as 

A, = [lnA,]/2i (4) 
A(-) ^ -[lnA(-)]/2t. (5) 

Obviously As is the negative Lyapunov exponent and xi 
is the positive Lyapunov exponent for forward time. It 
is easy to see that any theory for §00 and As would be 
applicable to s^'' and A'-"-', since they are defined for the 
same flow except for a reversal of the direction of time. 
From now on we concentrate on the set of (As, §00) and 
let 

X{tt) = -X,i^,t) (6) 
for convenience in notation. 

IV. THE HAMILTONIAN NATURE OF THE 
THEORY OF FINITE TIME LYAPUNOV 
EXPONENTS 

Before we proceed, it is useful to summarize the main 
results of ergodic theorem. Chaos means sensitive de- 



pendence on initial conditions. There are two aspects of 
chaos which are captured by the multiplicative ergodic 
theorem of Oseledec. The first one is the dynamical as- 
pect, i.e., how sensitive is the dependence on initial con- 
ditions? This is answered by the Lyapunov characteristic 
exponents. The second one is on geometry. Along differ- 
ent directions, neighboring points do behave differently, 
i.e., different Lyapunov exponent corresponds to differ- 
ent characteristic directions. As signified by its name, 
ergodic theorem treats the time-asymptotic limit or long 
time average. Hence the Lyapunov exponents in ergodic 
theorem are also called infinite time Lyapunov exponents 
which are constants in a chaotic sea. The geometrical in- 
formation that is given by the ergodic theorem is always 
local, i.e., §00 (0 ^ function of position. The smooth- 
ness of G(x) implies the smoothness of (C) as a vector 
field. The field line of ioo(0 of great importance in 
transport of advection-diffusion type. For example, the 
rapid diffusive relaxation of an externally imposed scalar 
or vector field occurs only along the s lines j?] . 

In applications where finite time is of concern, one 
needs to understand the properties of the finite time Lya- 
punov exponents. The finite time Lyapunov exponent 
has both a time and space dependence. It is also called 
local Lyapunov exponent for that reason. Since, for ex- 
ample, both X{^,t) and §00 (C) are local, one might feel 
that there could be some relationship which relates the 
spatial variation of these two. This was addressed by our 
work on finite time Lyapunov exponent 0|. We found 
that the finite time Lyapunov exponent A(^, i) can always 
be asymptotically decomposed into three main parts with 
the addition of an exponentially small fourth (correction) 
term 

A(e, t) = A(e)A + /($, t)/Vt + X^ + 0(cxp[-2A(^, t)t]), 

(7) 

where 

§oo-Vo/(e,t) = 0, (8) 
§00 • VoA(0 4- Vo • = 0, (9) 

and A°° is the infinite time Lyapunov exponent. We note 
that A(^) is a smooth function of position due to the 
smoothness of vector field Soo- 

The correction to the asymptotic decomposition, i.e., 
the fourth term in equation (^, becomes exponentially 
small as t becomes large. The rate of the exponential 
decay is given by the magnitude of the local Lyapunov 
exponent. Hence the correction term becomes negligible 
on a time scale of a few local Lyapunov time. It should 
be emphasized that what we presented here is an asymp- 
totic form of a local expression (finite time Lyapunov 
exponent is an explicit function of position and time). In 
different regions of a chaotic component, the local Lya- 
punov exponent can vary significantly. For example, the 
stochastic layer would have a much smaller finite time 
Lyapunov exponent and hence a longer period during 
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which the exponentially small correction is still impor- 
tant. However that time length is fixed on the order of 
a few local Lyapunov time, which is the characteristic 
time scale over which the stochasticity of a chaotic tra- 
jectory starting from a particular location is of practical 
importance. This should be contrasted with a statistical 
description of the finite time Lyapunov exponent, whose 
practical applicability would be affected by the extremely 
long transients for hamiltonian systems to reach to an 
invariant distribution by following one chaotic trajectory 

Because does not vary along the § direction, 

thus to exponential accuracy the variation of the finite 
time Lyapunov exponent along the stable manifold is de- 
termined by the geometry of the stable manifold alone, 
equation (O). Notice that everything is local, i.e. func- 
tions of position. In particular, A(^) and (^) do not 
have a time dependence. That is: they represent time 
asymptotic structures. They are of great practical im- 
portance because the equation accurately describes 
the spatial variation of the finite time Lyapunov expo- 
nent along the stable manifold on a time scale of a few 
local Lyapunov time. The effect of geometry on a dy- 
namical quantity like the finite time Lyapunov exponent 
is captured by the function A(^) alone, a function that is 
completely determined by the G(x) field. 

The main points of Tang and Boozer's work on finite 
time Lyapunov exponents are the function A and its 
relationship with §00. These results actually find their 
roots in hamiltonian dynamics. To see that, one can 
construct a new vector field 



S(?) 



.(0- 



Obviously S(^) does not vanish anywhere and is glob- 
ally defined in a chaotic region, or using Meiss's term, on 
an ergodic irregular component. More importantly, S is 
divergence- free because of equation . All the necessary 
information on chaos in field G(x) associated with nega- 
tive Lyapunov exponent for forward time, are contained 
in this new vector field S. 

On a regular component in phase space, i.e. KAM 
surface, S is also well defined and divergence-free. The 
property of a trajectory on a KAM surface is determined 
by the rotational transform L.^ip). In explicit form, the 
trajectory follows 

^ = ^0; = 00 + i'o(^)i; ^ 00 + L{i')Mi')t (10) 

with i'o{ip) the Jacobian of {ip,(j),9) coordinates. The 
form of equation ( p^ ) is generic for integrable trajectories 
psf . An explicit construction of the proper coordinate 
system which gives rise to equation (|l^) can be found in 
112[ . The §00 introduced earlier now takes the form 



(11) 



found in . Vanishing component of §00 means that 
Soo vector is tangent to the KAM surface. By definition, 
two neighboring points along §00 direction will converge, 
but quadratically in time, which should be contrasted to 
an exponential rate in a chaotic region. The S vector can 
be simply given by 



■i^o)- 



where prime denotes a derivative with respect to ip, the 
action. A detailed derivation of equation (11) can be 



This would correspond to an integrable trajectory with 
rotational transform {lvq)' /vq on the constant action ■(/^o 
surface. In other words, if one writes the original field as 

G = Vijj •kV0 + V(t) X Vx 
and on the constant surface 

i(-0o) = 9x/i9V'U=V'o> 
then the S field will be 

S = Vt/; X + V0 X Vx 
and the hamiltonian x for the field line of S satisfies 

5x/5V'U=Vo = (ti'o)7^'ol^=V'o- 

Cautions should be taken regarding the derivative with 
respect to ^. The surviving KAM tori are parameterized 
on a fractal set of action in a perturbed hamiltonian 
system. The proper definition of derivative with respect 
to V' on the surviving KAM tori invokes Whitney's no- 
tion, which was discussed in detail by Poschel Q . 

It should be noted that so far we have established 
divergence-free S field separately in chaotic region and 
KAM region. To achieve divergence-free globally, it ac- 
tually only requires the smoothness of the direction of S 
fields when crossing the last KAM surface. A jump in 
the magnitude of S is allowed since its gradient is per- 
pendicular to S (V • /Soo = V/ • Soo + /V - Soo). Although 
an analytical proof has not been found, numerical results 
have been obtained to support the continuity of the di- 
rection of §00 at the border between order and chaos. 

As pointed out in section ||, a globally divergence- 
free field admits a hamiltonian representation. Hence 
the field line of S (and s) is described by a hamilto- 
nian xi('0ij ^'ij </'i) with 01 the canonical position, ■01 
the canonical momentum, and 0i the canonical time. In 
other words, the characteristics of chaos in the hamilto- 
nian X are now contained in a new hamiltonian xi. Xi l^^s 
the same degrees of freedom as x and the field lines of S 
are chaotic just as those of G. Similarly, one could repre- 
sent chaos in S (xi) by another hamiltonian X2 associated 
with another divergence- free vector field S'. Henceforth 
a hierarchy of hamiltonians is constructed for describing 
the chaos in G. 

One might be concerned as to the proper counting of 
the degrees of freedom if G is time dependent. In actual- 
ity, S becomes time dependent through §00. This comes 
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from the fact that the specification of Lagrangian coordi- 
nates depends on the choice of initial time tg, as discussed 
in last section. For a time-dependent field G, the vector 
§00 is a function of position and time §00 (C; ^o)- If G(x, t) 
is periodic in time t, 800(^7 ^0) will be periodic in time to- 
The hamiltonian Xi has the exactly same degrees of free- 
dom as X- This small subtlety holds the key to a correct 
understanding of a time dependent field in two dimen- 
sions G{x,y,t) which has chaotic field lines. The time 
dependence in Sao{£,x, iy,to) assures that the correspond- 
ing S field has a time dependence and the hamiltonian 
Xi for S is of one and a half degrees of freedom. Hence 
chaos is allowed for the field line of Soc{^x,iy, h) and one 
s line fills the entire irregular component of G. 

As stated earlier, there is a symmetry between (§00, As) 
and (sio \ aI ^ ) . Hence one could construct a set of hamil- 
tonian to describe the chaos associated with the positive 
Lyapunov exponent and unstable direction for the for- 
ward time dynamics, in analogy to what we have done 
for the negative Lyapunov exponent and stable direction. 

V. SUMMARY 

In this paper we study the interrelations of divergence- 
free field, hamiltonian dynamics, multiplicative ergodic 
theorem in hamiltonian chaos, general coordinate trans- 
formation, and the theory of finite time Lyapunov ex- 
ponents. We argue that finite time Lyapunov exponents 
and its associated characteristic directions are the most 
important information on chaos for the purpose of phys- 
ical applications. Unlike the infinite time Lyapunov ex- 
ponents, spatial variation of the finite time Lyapunov 
exponent is directly related to the geometry of its cor- 
responding characteristic direction. Both are found to 
be given by one single new hamiltonian of same degrees 
of freedom. The magnitude of the phase flow field of 
the new hamiltonian determines the finite time Lyapunov 
exponent, while the direction of the phase flow gives the 
corresponding characteristic direction. The new hamilto- 
nian is chaotic on the irregular component of the original 
hamiltonian. We hope the point of view presented here 
could stimulate new insights into hamiltonian chaos. 
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